Phase Separation in Two-Dimensional Additive Mixtures 
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We study 2— dimensional binary mixtures of parallel squares as well as of disks. A recent cluster 
algorithm allows us to establish an entropic demixing transition between a homogeneously packed 
fluid phase and a demixed phase of a practically close-packed aggregate of large squares surrounded 
by a fluid of small squares. 
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Binary mixtures of impenetrable objects pose one of 
the important, and lively, problems of statistical physics. 
For many years it has been discussed whether objects of 
different types I (large) and s (small) would remain ho- 
mogeneously mixed as the number of these objects per 
unit volume increases. Particularly interesting cases con- 
cern socalled additive mixtures 0], as hard spheres with 
radii ri and or cubes with length di and dg. The prob- 
lem of phase separation in binary mixtures is of impor- 
tance as the simplest model for colloids. It has been a 
meeting ground for many different theoretical, compu- 
tational, and experimental approaches. As an example, 
the well-known closure approximations, as well as virial 
expansions, both of very great importance for the theory 
of simple liquids |Q, have been brought to bear on this 
problem, often with contradictory results. 



FIG. 1. Snapshot of 200 big squares {di = 1) and 20,000 
small squares {ds — 0.1) in a periodically continued box of 
size 30 X 30 (packing fractions rji = rja = 0.22). 

In this paper we discuss additive systems in two dimen- 



sions, parallel hard squares and also hard disks. For the 
parallel hard squares system, we clearly see a phase sepa- 
ration transition which has not been observed in previous 
work 1^,^. For hard disks, we expect an analogous tran- 
sition for extremely dissimilar sizes, which must satisfy 
rg/ri < 1/100 for packing fractions rji ^ rjs < 0.3. 

Monte Carlo simulations have long been performed on 
these systems. They were recently boosted by a new 
cluster algorithm [^|j6), which allows thermalization of 
systems orders of magnitude larger than previously possi- 
ble. The algorithm sidesteps a problem readily apparent 
in Fig. 1, which shows a typical configuration in the ho- 
mogeneous phase. There, each large square is surrounded 
by very many small objects. Trial moves of large squares 
will therefore be rejected in the overwhelming majority 
of cases, and the algorithm will get stuck quickly. Our al- 
gorithm rather swaps large patches of the configuration 
in a way which preserves detailed balance The 
algorithm is applicable for arbitrary shapes and in any 
dimension, and in continuous space as well as on the lat- 
tice. Most importantly, the method works even for ob- 
jects very dissimilar in size as long as the total density is 
not too high. 

We are able to converge our simulations of two- 
dimensional parallel squares for total packing fractions 
1] — rjs + r]i which do not sensibly exceed the percolating 
density rjperc- As in 3D, we notice that the density r]perc 
depends very little on the ratio R — ds/di < 1. We find 
rjperc — 0.5. Figs 1 and 2 show snapshots of the simula- 
tions for 200 large squares and 20, 000 small squares at 
equal composition [rjg — rji) and total packing fractions 
rj = 0.44 and rj = 0.60, respectively. It is evident that 
Fig. 1 represents a homogeneous mixture whereas the 
system shown in Fig. 2 consists of a 'solid' block of large 
squares surrounded by a fluid of small squares. In our 
opinion, these runs present direct evidence for a transi- 
tion of the homogeneously mixed fluid into a (solid-fluid) 
phase. 
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(periodically continued) lateral distances. With this def- 
inition, the distance of two large squares in contact is 
r ^ di and Gu = Vi i^^^(^\^\^\y\)^r9ii(x^y)dxdy with gu 
the usual pair correlation function In Fig. 4, one can 
see that the mixed system's Gn has pulled away from the 
pure system's correlation function on all scales, a clear 
indication that the large scale structure of the fluid has 
changed. The staircase structure is a very clear indica- 
tion of solid order. Finally, the reason for the rather large 
finite-size effects at large separation is immediately clear, 
since Gu{r) has to meet the curve of the corresponding 
monodisperse system with f] = rji for half the box length. 
A similar analysis was used to establish the instability 
line in Fig. 3. 



FIG. 2. Snapshot of 200 big squares (d; = 1) and 20, 000 
small squares (d2 = 0.1) in a periodically continued box of 
size 26 X 26 (packing fractions r]i — rja — 0.30). 

In the simulations at 77 = 0.6, we have slightly ex- 
ceeded the percolation threshold rjperc- Therefore, the 
algorithm will swap patches which usually comprise al- 
most the whole system. This generates problems for large 
systems, and we have e.g. been unable to converge (at 
r] = 0.6) a sample with Ni = 800, = 80,000. In con- 
trast, the simulations at lower packing fractions converge 
extremely rapidly for arbitrary system size. We can sum- 
marize the situation for equal composition (77^ — rji) by 
the diagram of Fig. 3: The gray area corresponds to the 
region of the diagram in which our algorithm performs 
extremely well. As mentioned, this region is delimited 
for the homogeneous system by the percolation threshold 
and by the appearance of high-density areas as a conse- 
quence of phase separation. 

We also studied the instability line as a function of 
the composition x = rji / iVs + Vi) j especially since re- 
cent work of Cuesta |^ suggests symmetry of the phase 
diagram with respect to x = 0.5 in the 3-dimensional 
system. In the 2-dimensional system, simulations at dif- 
ferent compositions {x = 0.3,0.5 and 0.7 for R ~ 0.1) 
contradict this behavior since critical packing fractions 
are 7?crit = 0.49 ± 0.02, 0.53 ± 0.02 and 0.60 ± 0.05, re- 
spectively. 

In our three-dimensional simulation , it was impos- 
sible to interpret the data by direct inspection, as in 
Fig. 1 and Fig. 2. We analyzed the transition therefore 
with the help of the integrated pair correlation functions 
Gu{r) — Airrji J r'dr'gu{r') We repeat the analysis in 
the 2D case in order to stress the soundness of our pro- 
cedure, which considers Gu rather than the much nois- 
ier gu. Gii{r) determines the average number of large 
particles around a given large particle within a distance 
r. For the special case of parallel hard squares, we de- 
fine r = max(Aa;, Ay), where Ax and Ay are the two 
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FIG. 3. Phases of the two-dimensional hard square sys- 
tem at equal composition r\s — rji (total packing fraction 
■q — Tjs + rji vs. R = ds/di). The squares locate the pa- 
rameters of the snapshots in Figs 1 and 2. The region in 
which our algorithm performs excellently is shaded in gray. 



A comparison of Fig. 4 with the data for the homoge- 
neous mixture (cf. Fig. 1) at = 0.44 is very revelling. 
In the latter case (not shown) , the difference between the 
mixed and the monodisperse system concerns mainly the 
region of small separation between the squares and would 
be unobservable on the scale of Fig. 4. The same obser- 
vation was made in 3D. The effects at small separations 
r are already detectable by visual inspection of Fig. 1, 
since there are quite many 'bound' pairs and triplets. We 
find that each large square has bound an average of 0.8 
squares. This agrees very nicely with our observation in 
3-D, where we noticed the onset of the phase-separation 
instability as the additional binding was close to one. 

We also performed simulations for mixtures of hard 
disks. Let us begin our discussion with an indirect heuris- 
tic argument which connects transition parameters for 
squares and for disks. It has long been understood 0| 
that the overlap of excluded volumes entropically favors 
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close contact of the large objects. The excluded vol- 
ume (for a small square) around a big square consists 
in the area of the latter and a strip of width ds/2 around 
it. Side-to-sidc contact between two squares leads to an 
overlap of excluded volume of the size AVgquare ~ x dg. 
As the large objects touch, the volume available to the 
small particles and therefore their entropy increase. At 
the same time, contact of large objects decreases their 
contribution to the entropy of the complete system. The 
phase separation transition appears when those two con- 
tributions are equals. As the decrease of entropy due to 
contact of large objects may be considered as indepen- 
dent of the ratio R of the size of the small and large 
objects, we may compare the increase of entropy due to 
larger available volume for small disks or squares. In fact, 
since the small particles' contribution to the entropy is 
S A^slogVs (with Vg the available volume), one finds 
A5 ~ (7V5/y)AV;quare, and therefore 
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(1) 



Repeating the same calculation for disks, we notice of 
course that the overlap of excluded volume AVdisk ~ 
\/rfri and the number of concerned small disks are much 
smaller. This leads to 
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FIG. 4. Integrated two-point correlation function for the 
dense system rji = rjs = 0.30 (cf. Fig. 2) for 50 and 200 large 
squares, respectively. The two curves are compared to the 
monodisperse system's correlation function. 

In conclusion we have studied the problem of phase 
separation of two-dimensional systems (hard squares, 
hard disks) by direct Monte Carlo simulation. For hard 
squares, our Monte Carlo data leave little room to doubt 
a direct fluid to (solid-fluid) transition. For hard disks, 
the stability of the homogeneous mixture seems estab- 
lished for any 'reasonable' ratio of radii r^/r; > 1/100. 
Our heuristic argument would however lead us to expect 
a transition for even more extreme ratios. 
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This order-of-magnitude argument tells us that (with a 
hard square transition at ds/di ~ 1/10 for i] ~ 0.5) we 
can expect a transition for disks at best for r^/rz ~ 1/100. 
Accordingly, our simulations for j-g/ri > 1/100 at 77 = 0.6 
have revelled no instability of the homogeneous phase. 
Even larger simulations, at r^/r; = 1/150, (50 large and 
1,250,000 small disks, 77 = 0.6) did not converge, even 
though the additional binding has continuously increased 
in the course of a month-long simulation. In these sim- 
ulations, the depletion potential is very strong, but also 
extremely short-ranged. The Monte Carlo simulation of 
such 'golf-course' potentials is of course extremely time- 
consuming and often impossible. 
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